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Abstract 


This paper discusses the rotating stall and surge in centrifugal compressors from the standpoint 
of vibrational diagnostics and rotor dynamics. The literature survey shows insufficient under- 
standing of rotor dynamic implications of rotating stall. Most of the literature on the subject 
addresses only aerodynamic aspects of the problem, thus putting an empress on stall cells with 
exclusion of the rotor vibrational response. This paper treats a compressor as a coupled 
mechanicakaerodynamic system. For the rotor/fluid system, the results of Dynamic Stiffness 
identification of an experimental centrifugal compressor in rotating stall and in normal operating 
conditions are presented. It is shown that the aerodynamically induced radial (direct spring) 
stiffness coefficients are negative in normal operating conditions. The negative radial stiffness 
effect in the case of rotating stall is even more profound. The result of it that the rotor system 
mechanical resonance could be reduced as much as to 5-10% of running speed. The conclusion 
is that the rotor dynamic stability margin during rotating stall is reduced. As a conclusion from 
the overview of case histories and aerodynamic researches, a diagnostic strategy is suggested. 


1. Introduction 


The interest in radial compressor fluid-induced instabilities started in the mid-seventies. Since 
then, a number of case histories of high-pressure centrifugal compressor instabilities have been 
published [1-17]. Most of the publications report two types of rotor vibrational behavior: 1) 
High eccentricity and rotor first natural frequency re-excitation, 2) Subsynchronous forward 
precession with rotative speed dependent frequency. The former is usually referred to as whip- 
type behavior and is normally associated with balance pistons, fluid film bearings, and labyrinth 
seals. The latter is called whirl type behavior [27-28] and can be associated either with fluid 
film bearings/seals or with rotating stall. Rotating stall usually results in an appearance of a low 
subsynchronous frequency component in the rotor vibration spectrum (frequency ratios fall 
typically between 8 and 40 %, but can go as high as 80% of rotative speed [20]). The emphasis 
on rotor behavior does not allow for details of particular flow patterns, but treats the fluid (gas) 
“in average,” which falls in the scale of rotor motion. 


2. Case History 
A 5-stage propylene compressor provides an example of the industrial centrifugal compressor 


behavior under the conditions of surge/rotating stall. The orbits (Fig. 1) show a typical forward 
precession with frequency 0.15x, supported by the lateral vibration, half spectrum cascade data 
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(Fig. 2, the data was taken before full spectrum becomes available). The gas pressure measured 
at the outlet of the third stage is presented as a timebase waveform for the steady state regime 
(Fig. 3) and spectrum cascade plot for the startup (Fig. 4). 
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Figure 2. Half spectrum cascade plot of the vertical 
lateral vibrations of the rotor at the outboard bearing of 
the 5-stage propylene compressor during startup. 


Figure 1. Typical orbits of the lateral 
vibrations of the rotor at the outboard 
bearing of a 5-stage propylene centrifugal 
compressor rotor in rotating stall 
conditions. 
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Figure 3. Gas pressure oscillations measured after the third stage of 5-stage propylene compressor in 
steady state regime. Note high vane passing frequency modulated by the low frequency of rotating 


stall/surge. 
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Figure 4. Half spectrum cascade plot of gas Figure 5. Half spectrum cascade plot of the rotor 


pressure oscillations at the outlet of the third 
stage of 5-stage propylene compressor during 
startup. 


axial vibrations of the Sstage propylene com- 
pressor during startup, showing surge component. 
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COMPRESSORS IN ROTATING STALL 
The steady state data shows the combination of vane passing high frequency with low frequency 
oscillations. Using the spectrum cascade plot, the latter are identified as surge component (4.2 
Hz) and rotating stall component (0.15x). It is important to note that the spectrum cascade of the 
thrust probe data (Fig. 5) shows a strong activity at the surge frequency, while the rotating stall 
component is much smaller. The described behavior is a reflection of the surge, which causes 
axial vibration at frequency 5 Hz at 2,000 rpm, while rotating stall is mainly causing forward 
circular radial vibration. 


3. Instrumentation and Diagnostics Strategy 


A typical machine train for a centrifugal compressor unit is shown in Fig. 6. It includes engine 
(in this case an aeroderivative gas turbine), gearbox, and the compressor itself. Different 
instrumentation setups for vibration, as well as performance diagnostics and monitoring of the 
engine, are beyond the scope of this paper. For the centrifugal compressor, the instrumentation 
package normally includes a Keyphasor® transducer, two sets of proximity probes in XY 
configuration installed at the inboard and outboard bearing locations (see Compressor DE and 
Compressor NDE in Fig. 6), two axial transducers, and performance monitoring instrumentation. 
As a minimum, the compressor performance monitoring instrumentation is supposed to include 
inlet and discharge pressure, temperature, and discharge flow measurements. Using the 
described instrumentation and case histories analysis, the rotating stall or surge can be identified 
based on the following information: 
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Figure 6. Typical compressor train diagram. 


e Radial vibration 
Steady state data: direct orbits, full spectrum 
Indication: subsynchronous forward component 
Possible causes: rotating stall, fluid- induced instability, and surge 
Transient data: direct orbits, full spectrum cascade 
Indication: subsynchronous forward component with frequency proportional to the 
rotative speed (0.1 to 0.8x) 
Possible causes: rotating stall, fluid- induced instability 
e Axial vibration 
Steady state data: spectrum 
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Indication: subsynchronous component 
Possible causes: rotating stall (possible weak coupling) and surge 
Transient data: spectrum cascade 
Indication: subsynchronous component with low constant frequency 
Possible causes: surge 
e Pressure, temperature, and flow 
Trend data: (dc signal) 
Indication: pressure head and flow decrease, temperature increases 
Possible causes: surge, less likely rotating stall 
Note: if pressure head and flow indicate a position on the performance curve close to the 
surge point, it could support rotating stall or surge versus fluid-induced instability. 
Dynamic data 
Indication: vane passing frequency amplitude modulation of pressure, with constant flow 
Possible cause: rotating stall 
Indication: pressure vane passing frequency amplitude modulation with correlated flow 
oscillations 
Possible cause: surge 


As indicated above, a strong conclusion can only be reached by the correlation of vibration data 
with performance information. It is important to note that rotating stall can originate in a 
particular stage of multistage compressor. To localize the rotating stall source, additional 
pressure transducers may be required. 


4. Experimental Compressor for Rotating Stall 


An experimental rig (Fig. 7) is constructed to evaluate Dynamic Stiffness for the rotor/fluid 
system during rotating stall and normal operating conditions. The rig consists of a single stage 
centrifugal compressor with a vaneless diffuser, driven by a l-hp dc electric motor. The motor 
speed is controlled by a high-output, variable power supply, and the speed can be controlled 
from near zero to 17,000 rpm. The motor is attached to the compressor rotor by a flexible disk 
coupling. 
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Figure 7. Experimental compressor setup. 
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The compressor rotor is separated into two pieces: a rotating section from the motor to the 
outboard end of the compressor and a nonrotating stub shaft from the outboard compressor 
bearing to the outermost support bearing. The rotating section is supported by two bearings: a 
pivoting bearing on driver-end and a rolling-element bearing on the nondriver-end in a soft 
support. A spring-load frame is attached to the nondriver-end of nonrotating stub shaft to control 
the radial position of the compressor rotor and counteract the gravity load of the rotor. The 
nonrotating stub shaft also provides the input for the perturbation set. The nonsynchronous 
perturbator transmits an exciting rotating force with variable frequency to the compressor rotor 
with a constant rotating force, 1.3 N. 


Two orthogonal pairs of eddy current displacement transducers are mounted next to the inboard 
bearing and at the spring-load frame to measure the displacement of the rotor. Two Keyphasor® 
probes are installed on the compressor shaft and the perturbator disk to provide phase and speed 
references. A conical exhaust valve on a lead screw controls airflow through the compressor. 
The valve is located at the end of the discharge pipe from the compressor (not shown in Fig. 7). 
It can be adjusted by several turns from fully closed, which completely blocks the output flow of 
the compressor. When the valve is fully opened, the airflow passes with minimum restriction. 
The actual flow is monitored by a flowmeter. Dynamic pressure transducers measure the 
compressor inlet and discharge pressure, as well as pressure wave propagation in the diffuser (at 
the same radius, 90° apart). 


4.1. Compressor performance 
An experimental compressor performance chart is displayed in Fig. 8 nondimensionally, 
applying general conversion rules as [27]. 


As it can be observed in Fig. 8 that with the airflow increase, the pressure drops, and the 
compressor is out of the rotating stall regime. Increase of the airflow stabilizes the compressor. 
The phenomena is supported by the full spectrum waterfall data, Fig. 9, taken at a constant 
rotative speed with the flow increasing by opening the exhaust value from full close to wide 
open. 
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4.2. Rotating stall experiments of compressor startups 

Under the flow constriction conditions, the experimental compressor rotor exhibits rotating stall 
whirl with a subsynchronous frequency component of 0.15x in the forward direction, presented 
in Figs. 10 to 12. Compared to the field case history above, the orbits (Fig. 10) show a clear 
forward precession with frequency 0.15x, supported by full spectrum cascade data (Fig. 11) 
generated from the lateral vibration. In addition, full spectrum cascade data (Fig. 12) is 
generated from pressure transducers, mounted in diffuser, further substantiating this result. 
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Figure 12. Full spectrum cascade plot generated from diffuser 
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operating conditions. 


4.3.  Nonsynchronous perturbation tests 

A series of nonsynchronous perturbation tests are performed at each of four compressor running 
speeds: 10,000, 12,000, 14,000, and 16,000 rpm. For the detailed description of the applied 
methodology see [29]. At each running speed, perturbation runs are performed under two 
different valve positions, chosen as close to each other as possible or more opened up for two 
values of mess flow. Those valve positions provide rotating stall operation or normal operation 
conditions respectively. The perturbator generates a force of magnitude F with orientation 5 
rotating in the same direction as the compressor rotation (forward perturbation), resulting in the 


perturbation vector Fe’. The perturbator is revved up from slow roll speed to 5,000 rpm in 
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order to identify Dynamic Stiffness of the rotor first mode. The vibration signals from two 
orthogonal, eddy current displacement transducers located at spring-load frame as 45° (x) and 
135° (y) to the left (as seen from the driver) are filtered to the frequency œ of the perturbation 
force. This yields amplitudes and phases for vertical (y) and horizontal œ) components of the 
response. 


In order to minimize the influence of mechanical stiffness asymmetry, 1x vertical Y, and 


horizontal X, vectors are combined into the circular forward (in the direction of perturbation) 
component: 


AM) =¥, + jXq (1) 


Next, the perturbation process is repeated for the perturbation in a direction opposite the rotation 
of the compressor rotor (reverse perturbation) for each running speed Q and valve settings. In 
this case, the forward (in the direction of perturbation) component of the filtered-to-perturbation- 
frequency response is opposite the direction of rotation. The nonsynchronous Dynamic Stiffness 
is calculated by taking the ratio of the force input to the response vectors for each speed in the 
databases. 

Fe? 
DS (w)= DDS(w)+ jODS (w)= FUR 


@ 





, DDS (o)=direct(DS), ODS(w)= quad(DS) (2) 


Figures 13 and 14 present the forward and reverse perturbation data with the results of Dynamic 
Stiffness under the test conditions described above. The negative speed there corresponds to the 
reverse perturbation. 
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Both figures show that in the forward direction, the direct Dynamic Stiffness increases with the 
rotor speed increases, but in the reverse direction, the direct Dynamic Stiffness keeps about the 
same values for all rotor speeds. This phenomenon indicates the effect of the fluid inertia. 


Since the total Dynamic Stiffness of the compressor consists of mechanical and aerodynamic 
parts, the former is measured by running nonsynchronous perturbation on the compressor at zero 
rotative speed (DSo). The difference of Dynamic Stiffness between the running compressor and 
stopped compressor represents aerodynamically induced Dynamic Stiffness, DSair: 

DDS ,,,(w) = direct(DS — DS, )= K.,, —M ,,,,.®° —M,,,,.(@ -AQF (3) 
ODS ,,, (o) = quad(DS — DS , ) = D,,,(@ —2Q) 
where Mairi , Kair, Dair are air inertia, radial stiffness, and damping respectively, generated by the 
air/impeller interaction. Mai is a fluidic mass effect term, and À represents a swirling ratio 
coefficient. 


Equation (3) shows the parameters of direct and quadrature Dynamic Stiffness extrapolation. 
The results of the raw data and the described processed data are illustrated in Figs. 15 and 16. 
Figure 15 shows the reductions of air-induced radial stiffness. The swirling ratio coefficient A is 
shown in Fig. 16. The negative value of À indicates the reverse-flow pattern rotation. It gets 
stronger at lower rotor rotative speed and in the rotating stall operating conditions. 
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Figure 17. Eigenvalues for rotor in rotating stall and normal 
conditions comparing to mechanical model. 
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Equation (3) allows also performing eigenvalue calculations for the compressor rotor dynamic 
stability evaluation, shown in Fig. 17. The eigenvalue analysis clearly shows that the rotating 
stall lowers radial stiffness and decay rates. In the case of rotating stall, increase of the rotative 
speed increases the rotor stability margin. The result of eigenvalue analysis shows that rotating 
stall not only lowers the rotor radial stiffness, but also brings the compressor toward the edge of 
the instability. (In the Fig. 17, the negative frequencies indicate the reverse perturbation). 


Another set of perturbation tests is performed for increased the airflow of the compressor. The 
tests also include forward and reverse perturbations with two rotative speeds of the compressor 
rotor. In order to produce the Dynamic Stiffness (Fig. 18-B) from the perturbation data (Fig. 18- 
A), the same routine is used as Eqs. (1) and (2). The same concept of subtracting the mechanical 
Dynamic Stiffness from total Dynamic Stiffness (Eq. (3)) is used to determine the aerodynamic 
part of Dynamic Stiffness. The results are illustrated in Figs. 19 and 20. They show that the 
airflow increases the direct Dynamic Stiffness (Fig. 18). This offsets the negative effect of the 
radial stiffness induced aerodynamically (Fig. 19). The swirling ratio coefficient à becomes 
positive as the airflow increases (Fig. 20). Figure 21 shows the result of the corresponding 
eigenvalue calculations, showing that when the airflow gets stronger, the compressor moves 
away from rotating stall regime and becomes more stable. 
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5. Conclusions 


This paper provides some practical recommendations on the rotating stall diagnostics as well as 
some insight into rotor dynamic implications of the phenomenon. It fills the gap in the literature 
on rotating stall, which all but ignore the rotor dynamic implication of rotating stall. The most 
significant conclusion from this paper is that aerodynamic forces applied to the impeller generate 
negative radial stiffness (stiffness Kair appears negative in all experiments). It is also important 
that this stiffness drops even more during the rotating stall. In short, the conclusions can be 
formulated as follows: 


e The rotating stall conditions manifest in the rotor hteral vibration signature as rotor forward 
precession with subsynchronous frequency, which tracks the rotor speed. This behavior is 
referred as a rotating stall whirl. 

e From the diagnostic standpoint, rotating stall differs from other instabilities of the whirl type 
in its strong dependence on the compressor operating conditions; it should disappear with 
increased flow. It differs from surge by its frequency proportionality to the rotative speed as 
well as by the predominantly axial character of surge-caused vibrations. 

e Surge conditions manifest in the rotor response as predominantly axial vibrations with low 
frequency independent of rotative speed. 

e The parameter clearly identified by nonsynchronous perturbation tests is aerodynamically 
induced radial stiffness. In all conditions, it is negative. From the practical standpoint, it 
means that even in normal operation conditions, centrifugal compressors have lower natural 
frequencies than it appears in purely mechanical calculations. 

e The aerodynamically induced radial stiffness during rotating stall drops dramatically and 
ultimately can make a compressor unstable. The data from field case histories indicates that 
the resulting mechanical resonance can be driven by rotating stall as low as 5-10% of running 
speed. 

e Altogether, aerodynamic forces may have a significant destabilizing effect on the 
compressor. In the case of a flexible rotor and/or soft bearings, they are able to create a rotor 
dynamic instability of whirl type. 
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